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1. Introduction 

Heavy ion reactions at intermediate energy offer great promise to measure spectro- 
scopic properties of nuclei far from stability, but one needs a tractable reaction theory 
to interpret the experiments. In this respective the availability of higher energy heavy 
ion beams is most welcome, because it becomes a reasonable theoretical approximation 
to neglect exchange of nucleons between the colliding nuclei. One can therefore consider 
the interaction in each nucleus as that of an external (complex) potential field. Within 
the framework of this basic approximation and with given potentials, we will here address 
the question of the accuracy of further simplified models of the reaction cross sections. A 
number of theoretical models have been proposed and calculated []I|-|T2[ in which different 



approximations were made. In this work we focus on two of the models, the eikonal model 
[0]"[0 transfer-to-the-continuum (TC) model 1§]-[|T3 . 



There are several cross sections that are measured and calculated in the models. The 
simplest measurement is the single-neutron removal cross section, in which only the pro- 
jectile residue, namely the core with one less nucleon, is observed in the final state. Besides 
the integrated removal cross section, denoted by a-n, the differential momentum distri- 
bution (Pa-n/dk^ is also measured. A particularly useful cross section is da-n/dkz, the 
removal cross section differential in longitudinal momentum. If the final state neutron can 
also be measured, the corresponding coincident cross section Ap {Ap — 1) + n is called 
the diffractive breakup cross section. The difference between the removal and diffractive 
breakup is called the stripping cross section. 

This paper is organized as follows. In Sec. II we summarize the essential ingredients 
of the TC model and the eikonal model, and discuss the accuracy of the neglect of finite 
interaction times in the eikonal model. In Sect III we discuss the different treatment of 
the neutron-target S-matrix in the two models. We will not discuss the accuracy of the 
Hankel function approximation separately, but it of course plays a role in the comparison 
of cross sections that we make in section IV which contains also our conclusions. 



2. Theoretical models 

All theoretical methods used so far rely on a basic approximation to describe the colli- 
sion with only the three-body variables of nucleon coordinate, projectile coordinate, and 
target coordinate. Thus the dynamics is controlled by the three potentials describing 
nucleon-core, nucleon-target, and core-target interactions. In most cases the projectile- 
target relative motion is treated semiclassically by using a trajectory of the center of the 
projectile relative to the center of the target s(t) = bc + vt with constant velocity v in the 
z direction and impact parameter in the xy plane. Along this trajectory the amplitude 
for a transition from a nucleon state ipi to a state ipf is given by 

1 

Af^ = — dt< ijf{t)\Vnt{v)\ij,{t) >, (1) 

III J— CO 

where Vnt is the neutron-target interaction. The state ipi will be the bound state of the 
nucleon in the projectile, while the final state ipj is a. continuum state. The detailed 
derivation of Eq.(|l|) from a scattering amplitude containing the full time dependent prop- 
agator can be found in Sec. II of Ref. |^. There it was shown to hold under the hypothesis 
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that the breakup process is hmited to peripheral projectile-target trajectories and that it 
is due mainly to the neutron interaction with the target potential. The probabilities for 
different processes can be represented in terms of the amplitude as 

^ = EI^^Ne-e/), (2) 

where ^ can be momentum, energy or any other variable for which one measures a differ- 
ential cross section. 

The effects associated with the core-target interaction will be included by multiplying 
the above probability by the 6c-dependent probability for the core to be left in its ground 
state. 

Thus the differential cross sections with respect to the longitudinal momentum is 
^^C^sf,^/-^^P.M. (3) 



where kz is the (longitudinal) recoil momentum of the neutron (see Eq. (2.3) of |TT|) and 
C^S is the spectroscopic factor for the initial single particle orbital. The cross section 
can be further divided into a stripping cross section agtr and a diffractive breakup cross 
section Cdiff depending on whether the removed neutron is detected in the final state or 
not. We accordingly will consider these individual probability distributions, 

P-n = Pstr + Pdiff, (4) 

and use a similar notation for the cross sections. 

We first summarize the transfer-to-the-continuum (TC) model. This model treats the 
time-dependence of the reaction explicitly, thus conserving energy. It uses the on-shell 
neutron-target scattering matrix, therefore making it in principle model-independent with 
respect to that interaction. It also makes use of the asymptotic form of the neutron wave 
function in the projectile. This is an asset in that the formulas have an analytic limit, 
but a disadvantage in that the results are only reliable at peripheral impact parameters]^. 
The two breakup probabilities are given by the following expressions: 

dP^Ufibc)_J2\^_S^jB{L,kzX), (5) 



dkz 

dPstrjb, 
dk. 



^(l-|^,j2)i?(/„,fc„6,), (6) 



where the factor -B(/n, ^2, &c) is a transfer probability which depends on the details of the 
initial and final states, and on the energy of relative motion. It is given by 

S('».t.M = 2(-)^(2'..+ l)|C.r^M,,,, (7) 



^ The derivation of the TC model requires that the neutron-target and the neutron-projectile potentials 
do not overlap. In this respect the result that can be expressed entirely in terms of asymptotic properties 
reminds one of Beg's theorem []l3|. 
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where k is the angular momentum of the bound neutron in the initial state with respect to 
the core. The variable has the interpretation as the angular momentum of the neutron 
with respect to the target. Also 



M, 



1 7 

/ dxe-^ PiXX, + B,x^)PiSXf + BfX^), (8) 
a/tt jo 

The arguments of the Legendre polynomials Pi. and Pi^ are Xj = 1 + 2{kz/'^iY and 
Xf = 2{{kz — mvY/{^f + 2mvkz — {mvY) — l , Bi = 2i]/d'yi and Bf = 2r]/dkf. 7j is related 



to the initial state binding energy by 7^ = \/ —2mei/fi. The variable V = y^z ~^ 1i has the 
interpretation as the modulus of the transverse component of the neutron momentum. 
In our notation is the energy of the neutron relative to the target in the final state. 
For diffraction this is the same as the final laboratory energy of the neutron if the target 
recoil kinetic energy is neglected. In the case of stripping, if it goes through compound 
nucleus formation, is the excitation energy of the compound state above the neutron 
threshold in the residual nucleus. For inelastic scattering it is the energy of the breakup 
neutron before it scatters from the target. This is equivalent to the sum of the excitation 
energy of the target final state and the final neutron energy relative to the target. If the 
target recoil kinetic energy is neglected the final kinetic energy Ej of the ejectile is given 
by the energy conservation condition 

Ef — Einc = Q = Si — Ef, (9) 

where Ei^c is the initial incident energy of the projectile in the laboratory, Q is the 
reaction Q value given hy Q = Si — e j and Si is the initial neutron binding energy in 
the projectile. With this approximation Eq. (|^) relates to the the projectile residue 
parallel momentum. Finally, |Cjp is the asymptotic normalization constant of the initial 
bound wave function 

^,{r) = -^'Ca^hf^{^l,r)Yl^^XOA), l^r » I. (10) 

It is obtained by fitting a realistic radial wave function to the Hankel form h\^^ outside the 
potential radius. In this way the transfer to the continuum results are model dependent. 
On the other hand the dimensionless quantity A/^ = 7~^|CipC^S' has been called reduced 
normalization in the contest of spectroscopy done with transfer reactions [ P^p3[ . If one 
considers A;, as the prefactor of the theoretical cross section obtained with the Hankel 
function, the ratio between the experimental cross section and the theoretical cross section 
would determine its value. 

Eqs.(|^-|^) were derived assuming no overlap between neutron-core and neutron-target 
potential. This assumption can be avoided if one makes an eikonal approximation to the 
basic expression for the amplitude, Eq.(|l]), as shown in Appendix B. Then one derives 

1,0 



dkz 27r JO 

where bn is the transverse coordinate of the neutron with respect to the target. The 
neutron-target S'-matrix is approximated by the eikonal form S{hn) = e"*-*^*-^"-*, related to 
the optical potential Kit by 
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X(bn) = ^ j^^Vnt{x,y,z')dz' . Finally, \i)i{hn -hc,k^)\^ is the longitudinal Fourier 
transform of the initial state wave function. Also the connection between Eqs.(^0) and 
Eq. (0) is made by replacing the sum over partial waves in Eqs. by the integral over 



impact parameters as in Eq. (|TT]) , and by evaluating the longitudinal Fourier transform in 
Eq.(0) using the asymptotic initial state wave function, Eq . (|10D . The one-dimensional 
Fourier transform of the initial wave function can be calculated analytically in the case 
of an Hankel function approximation Eq. ([lO|) yielding 

[Zli + i) ^. [Zli + i) 

^ Cf^PiXX^), (12) 

27]p 

where p = |bc — bn|. 

The total breakup probability is obtained from the integral of Eq.(|lTl) involving 



f^mm ^j^^ f^max kiucmatically allowed minimum and maximum neutron parallel mo- 

menta discussed in the following. Therefore although Eq. (^) describes the neutron-target 
rescattering in the eikonal approximation, it still satisfies neutron energy and momentum 
conservation. 

On the other hand, a sudden approximation or frozen halo form of the eikonal method 
has been used in where energy conservation was neglected. This approximation 

can be derived from Eq.(ll) when the integration limits in Eq.(^) can be extended to 
±oo. Then Eq.(0) is just be the longitudinal density, and the formulae for the TC and 
eikonal model become identical. In fact in this limit the removal cross section reduces to 

a.n = C'S J d'Kj dh^ [1(1 - S)\' + 1 - \S\'] |5,t(&c)lV.(bn - be,^)P, (14) 



which is consistent with the breakup cross section of In Eq.(^4D 5* is an eikonal 

S-matrix, as defined after Eq. ([Tl|) . 

The steps necessary to obtain Eq.(0) from Eq.(|ll|) and Eqs.(|-|^) can be justified in 
the high energy limit as follows. If the neutron binding energy is not too large, the 
final energy or momentum distributions are strongly peaked at the incident energy per 
nucleon. Therefore it is possible to average the dependence of the neutron target optical 
model S-matrix of Eqs.(^,^ over the full range of neutron continuum energies and assume 



that the S-matrix can be approximated with the eikonal values S obtained at Ef = |mu^ 
On the other hand, for the eikonal approximation to be good it is necessary that the 
parallel component of the neutron momentum kz be large with respect to the transverse 
component rj. In particular such condition must be satisfied for the minimum values 
of both. From the definition of k^ it is easy to see that its lowest possible value is 
j^mm _ ^mv'^)/{hv) in correspondence to = 0. And that the minimum value 

of ?7 = 7j. Therefore we get >> 7^ if a/Ej^c/S >> ^/ei. For a real halo with 

separation energy around 0.5MeV the condition >> 7j is satisfied at all initial 
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energies. Increasing the binding energy it is necessary to go to higher incident energies 
in order for the parallel momentum component to be larger than the transverse one. For 
a typical binding of lOMeV or more the conditions for the sudden eikonal approximation 
to be valid are satisfied from about 80A MeV. It is useful also to consider the values 
of neutron-target center of mass momentum fc^ = ^v/h where /i = At /{At + 1) is the 
neutron-target reduced mass. Here too we see that for a typical diffuseness of a=0.5fm 
the semiclassical condition ak^ >> 1 is satisfied starting from about 80MeV. Finally if 
such condition is satisfied we can extend the lower limit in the integral to —oo. There 
is however also an upper limit to the final neutron or ejectile momentum value in the 
parallel direction, due to the energy and momentum conservation. It was discussed in 
ref.[jT6[ where it was shown that 

m is the neutron mass, k^""^ = 0.96/m^^for a ^"^Be projectile. In terms of the maximum 
neutron final energy in the continuum the above condition reads 

^ 2 

If this upper limit is also extended to +00 then the kz integral gives a 5-function and 
inserting Eq. (|TT|) and the damping factor = |>S'cjp in Eq.(|I) we finally get Eq. ([I^) . 

These very strict conditions necessary to extend the theoretical limits of the k^ integral 
to infinity, can be somewhat revised in practical calculations. To see how accurate the 
sudden approximation is, we have calculated the integral J(/c^™, fc™''^) under various 
conditions of angular momentum, neutron binding energy in the projectile, and projectile 
velocity. For values of the parameters of interest there can be a rather large reduction 
for small values of the neutron transverse radius in the projectile, \bn — bd- However, the 
approximation becomes increasingly accurate as the transverse radius is made larger. In 
Fig. 1 we show the regions of energies in which the ratio 

^ /(-oo,oo) 



is within 5%, 5-10% and 10% of being unity, for a ci-wave orbital at a transverse radius 
corresponding to strong absorption. The region of good agreement is somewhat larger 
for s- and p-waves. Thus it seems that the sudden approximation should be acceptable 
for weakly bound neutron orbitals (< 6MeV) and beam energies greater than 40A MeV, 
provided the reaction samples the initial wave function at distances rather larger than the 
sum of the radii of the reaction partners. 

It must be noted that here as in Eq.@ we have adopted the so called no-recoil ap- 
proximation 1^ which consists in factorizing out from the matrix element the core-target 
S-matrix. Another important remark has to do with the extension of the lower limit in 
the kz integral to —00. Values of kz smaller than fc™" obtained when Sf = would 
correspond to final states in the projectile below the breakup threshold for the diffraction 
term. For the stripping term they would correspond to final bound states in the target. 
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Taking fc™" — oo is analogous to the completeness relation introduced in For this 
reason in Eg. ([T^) the diffraction term has usually been corrected by subtracting terms 
like I < — SSct\'4^i > In the transfer to the continuum model this correction is not 
as important because the neutron final energy is always positive and the integration is 
performed in the kinematically allowed region. 

Estimates of absolute breakup cross sections published so far have been made either 
using Eq.(^ with the breakup probability given by the transfer to the continuum the- 
ory Eqs.(|-H) IIT0H12I or with Eq.(0) |jl]^J^J^,|l^ . There are a number of assumptions 
contained in both models, the vahdity of which we are going to study in this paper. 



3. n-target optical potential and S-matrix 

Experimentalists often use light targets such as ^Be or ^^C to study spectroscopy by 
projectile fragmentation reactions. The definition of an optical potential for light targets 



is a very delicate issue which has been discussed in the literature for long time |]18[. There 
are two main issues, i) global parameterizations are based on large nuclei; ii) the imaginary 
optical potential changes drastically for light nuclei. Recently these problems have been 
faced in the contest of halo breakup studies ||5| JTl|Jl^ where different choices have been 



made by different authors. Bertsch et al.[^ used the Varner et al. parameterization 



rof . Such parameterization has been obtained taking into account relatively low energy 
proton and neutron cross sections (10-26MeV) on heavy targets, A=40-209. In ref. |^ 
a neutron target S-matrix was constructed in the optical limit of Glauber theory. The 
corresponding optical potential does not need to be given explicitly, but the behavior of the 
transmission coefficient given in would suggest a strongly absorbing potential of volume 
type. Bonaccorso fill] extended to high energies a parameterization of a phenomenological 
n-^Be optical potential fitted to low energy data. This potential seemed however to 



overestimate the free particle cross sections at high energies |]TT| and therefore ref.|T2 
also considered a microscopic optical potential calculated according to the method of 
Jeukenne, Lejeune and Mahaux (JLM) [jl8|. The JLM potentials are more complicated 
to calculate and to use than a simple standard parameterization. Then for the purpose of 
the present paper we have attempted to fit the available experimental n-^Be total cross 
sections with a new potential of the standard form. The potential is given by a real 
Woods-Saxon well and it has both volume and a surface derivative Woods-Saxon forms 
for the imaginary part. In obtaining the parameters given in Table I we have been guided 



by some existing parameterizations |2O|,0 obtained by fitting low energy data. We have 
modified such parameterizations to get a smooth behavior of the free particle neutron-^i?e 
cross section in the energy range 10-180 MeV. The free-particle cross section have been 
obtained in an optical model calculation. A good agreement of calculated free neutron 



angular distributions with the data of ref. has also been obtained. This potential 
satisfies very well the subtracted dispersion relations as given for example in Ref.[^ at 



energies larger than 40MeV. For lower energies the agreement is less good at small radii. 



In Fig. 2, top left, we show the experimental cross sections p3| together with the optical 
model cross sections, where the full curve is the total cross section, the dot-dashed curve is 
the elastic cross section while the dashed curve is the reaction cross section. At the bottom 
left we show the same quantities calculated in the eikonal approximation. The effect of the 
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eikonal approximation is to reduce the calculated cross sections, in particular at the lower 
energies. The reduction is more pronounced for the elastic cross section, probably because 
quantum mechanical reflection effects on the potential barrier can not be reproduced by 
the eikonal approximation. For comparison we show at the top right and bottom right 
respectively the same quantities calculated with the Varner parameterization. Varner 
potential does not reproduce the high energy data, also it gives a dominance of the elastic 
scattering up to about SOMeV. Here too the eikonal calculations underestimate the elastic 
scattering at low energies. 

Our potential has a rather strong surface term. It gives also large volume integrals 
{Jw/^ ~ 200MeV"/m^) in accordance with the light nuclei systematics [|T^ , P^] . The mi- 



croscopic origin of it can be understood since '^Be is weakly bound and it has rather high 
breakup probability. Also it has been known for long time pSf to have an unusually large 
mean square radius and to be strongly deformed. Finally we notice that we get a very 
large elastic vs reaction cross section ratio at low energy. The JLM potential gives the 
same behavior. The little experimental free particle data available show similar trend and 
it would be interesting to see if breakup reactions around 20A MeV which will soon be 
feasible at GANIL, will parallel such a behavior. 

In Fig. 3 we show the behavior of the term |1 — 5*1^ and of the transmission coefficient 
1 — when calculated by the optical model (solid and dashed line respectively) and 
by the eikonal approximation (dotted and dot-dashed line). The results are shown as 
a function of the angular momentum /„ and in the case of the eikonal calculations the 
semi classical relation /„ -|- 1/2 = khn was used to make the connection between the angular 
momentum and the neutron impact parameter. We give results at energies ranging from 
20 to 100 MeV. An interesting characteristic of these results is that the eikonal seems to 
concentrate the scattering at larger impact parameters than in the optical model case. 
Such effects have been seen by other authors ||2^ in different situations that ours and 
sometimes it has also been shown [0,^ that higher order eikonal corrections can improve 
the agreement with the quantum mechanical calculations. It is possible that in our case 
the effect is amplified by the deep surface part of the imaginary potential. It is well 
known that the eikonal approximation works well in the presence of strongly absorbing 
volume potentials and differences would also be less extreme for potentials that have more 
diffuse edges. This is because the reflection effects at the barrier that modify the neutron 
trajectory are not taken into account in a semiclassical approach. 

From the differences seen here, we would expect quite large differences between predic- 
tions of the two treatments of the S-matrix. The quantiflcation of these differences is the 
subject of section IV. 



4. Cross sections and conclusions 

We start this section by showing some results for the total probability 
Pipe) = P^n{pc)Pct{pc) obtained from the integrand of Eq.(|^) after integrating over kz and 
using Eq.(|18]). The Rg values are given in Table III at each energy. They were obtained, 
as explained in Appendix A, from the S-matrix calculated by folding the neutron-target 
optical potential of Table I with the projectile density, at all energies but 20A MeV. At 
such a low energy the optical limit of the eikonal model cannot be justifled, but calculations 
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by Carstoiu in second order eikonal approximation, folding the JLM potentials agree 
with our Rs at all energies, including 20A MeV. Same values are also obtained from the 



optical potential of |28 



For our numerical comparison of the various methods, we have chosen to study the 
breakup from the s and p states in ^"^Be which were recently measured by Navin et 
al.|T^. The initial bound state parameters are given in Table II. 

In Fig. 4 we show the total probabilities for diffraction obtained in the TC calculation 
by the solid lines. The dashed lines are for the absorption or stripping. Energies are as 
in Fig. 3. The dotted lines are the calculations from the diffraction integrand of Eq.(|lT]) 
while the dot-dashed lines are for the stripping. Because of the shift in the S-matrices 
of Fig. 3, we see here that the eikonal calculations have their maxima at slightly larger 
core-target distances than the optical model calculations. At the lowest energy the optical 
model gives almost equal diffraction and stripping probabilities. The eikonal calculation 
gives a rather larger diffraction probability. In all other cases the stripping is the dominant 
term. This effect was first pointed out in normally bound nuclei P^,P3 and it has been 



confirmed by some experiments and other theoretical calculations [[12 |. 

The TC calculations discussed in this paper are performed using the asymptotic, Hankel 
form of the initial state wave function, which can be the origin of some divergence if 
the projectile and target potential have an important overlap during the reaction. This 
problem was discussed in great detail in Ref . where it was shown to be more important 
for the stripping term since it is proportional to the neutron target phase shift, while 
the diffraction term, being proportional to the square of the phase shift converges more 
rapidly. This situation reflects the fact that diffraction reactions involve a more peripheral 
part of the initial wave function than stripping reactions. Here we can easily see from 
Eq.(0) that the integral over 6„, the neutron target impact parameter is well behaved 
provided the neutron-target potential and the corresponding S-matrix fall off more rapidly 
than the tail of the initial wave function. A practical way to satisfy such a condition and 
to avoid any divergence in the sum over neutron partial waves in Eq. (^) is to use a small 
value for the diffuseness of the imaginary part of the optical model (cf. Sec. IV of 0). 
Otherwise, as it is shown in Appendix B it is possible to start again from Eq.(|l]) and 
obtain a formula which contains the realistic single particle wave function and an optical 
model S-matrix, still keeping the information on the proper kinematical limits for the 
neutron parallel momentum. The only hypothesis involved is that high neutron angular 
momenta are dominant. 

To show the amount of error that is introduced by the Hankel approximation in the 
TC method, we show in Fig. (5) the integrand functions of the stripping and diffraction 
terms of Eq. (|lll) , after integrating over kz, obtained from the realistic bound state wave 
function and the corresponding terms in the sum over partial waves of Eqs.(|^J^) in the 
case of the incident energy of 78A MeV. The calculations were done at the fixed impact 
parameter = b.Qfm between the projectile and target. Left figure is for the diffraction 
term while right figure is for the stripping term. We can see that the calculation for 
stripping according to Eq. which is done with the Hankel function and optical model 
S-matrix (diamonds) shows a slower decrease in the interior of the projectile ( small 
hn ) than the calculation with the realistic wave function (crosses) done according to 
Eq.(^) of Appendix B. Therefore the comparison between these two results gives an 
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indication of the effect of tfie Hankel function approximation. On the other hand the 
eikonal calculations (full curve), done with the full wave function have higher maxima 
which occur at a larger distance than the optical model calculation. This effect is due 
to the behavior at large impact parameters of the eikonal S-matrix, as we have discussed 
before. The comparison of the full curve with the crosses gives an indication of the effect 
of the eikonal approximation. In the end the integrated probabilities obtained with the 
three possible approximations contained in Eqs. and (H), Eq.(^) and Eq.(0) differ by 
not more than about 10%. This leads to about 20% difference in the total cross sections 
because of the further integration over be, the core-target impact parameter. From the 
peak values of the diffraction and stripping curves it is easy to see that the stripping is 
almost three times larger than the diffraction. 

The values of the total cross sections are given in Table III. Because the eikonal proba- 
bilities are shifted towards larger distances we get total cross sections which are larger by 
about 20% than the TC calculations at all energies but 20A MeV. At this lowest energy 
two effects tend to compensate each other: the eikonal free particle cross section is largely 
underestimated while there is an overestimate of the phase space kinematically allowed. 
Then the eikonal and optical model calculations seem to give very close results. 

Considering all the effects discussed above we conclude than present models used to 
calculate the breakup cross section might tend to overestimate the true value, but the error 
should be of the order of 20% or less. On the other hand presently available experimental 
data are inclusive with respect to the target which was unobserved. In the case of ^Be it is 
possible that some part of the cross section comes from reactions in which the target itself 
underwent breakup. These reactions could account for about 20% of the measured cross 
sections and they are certainly not taken into account by presently available theoretical 
models. We suggest therefore that future experiments detect the final excitation state 
of the target in coincidence with breakup events (as done for example in Ref. |^ for 
"normal" heavy ions) and at the same time that theoretical models be improved through 
a more realistic treatment of the core-target interaction, as mentioned in Sec. IV and in 
Appendix A. We propose also an improved calculation of the breakup amplitude using 
the new formula Eq.(^) derived in Appendix B that combines the better treatment of 
the wave function in the eikonal model with the better treatment of the target interaction 
in the transfer-to-continuum model. We leave to a future study to apply this equation to 
the reaction data. 

Ejectile momentum distributions are not shown in this paper. They have been discussed 
in great detail in |10-1^ where it was shown that there can be noticeable differences in 



the results of the TC and eikonal methods. Present data do not have enough statistics to 
distinguish clearly between the two models. 
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A. Ion-ion S-matrix 

We discuss here the theory of the core-target S-matrix, needed to calculate Pctipc) in 
Eq.(^). The best source of information comes from the ion-ion reaction cross section 
which is related to the core-target S-matrix by 

at = (fh,[l-\S,tih,)\'). (17) 

Here we have assumed that the semiclassical replacement of the partial wave sum by 
an integral over be is permitted. Unlike the nucleon-nucleus scattering, this is always 
a good approximation, because of the smaller wave length involved, and one can safely 
calculate this S-matrix in the eikonal approximation. The usual procedure is to define 
an optical potential for the core-target scattering, and calculate the S-matrix from it. In 
practice the behavior of the S-matrix is determined by two parameters. The first and 
most important is the strong absorption radius Rs, defined as the distance of closest 
approach for a trajectory that is 50% absorbed from the elastic channel. The reaction 
cross section practically speaking is determined by the strong absorption radius. The next 
most important parameter is the thickness a of the absorption region. In optical model 
fits, it is closely related to the asymptotic behavior of the imaginary potential W{r) ~ 
gC-?"/"), Just as tJie reaction cross section strongly constrains the strong absorption radius, 
measurements of the elastic angular distribution provide information about a. 

In this work we shall use a simple parameterized form for Pct(bc), namely 

Pet(&c) = exp(-ln2e(^=-^^)/''). (18) 

It would be correct to use the elastic scattering S-matrix if the experimental cross 
sections were reported for projectile breakup leaving the target in its ground state. In 
fact the target final state is difficult to measure and most experimental cross sections are 
inclusive with respect to the target final state. Thus the measured cross sections should 
be somewhat larger than one would calculate with an elastic core-target S-matrix. 

A way to include absorption of the projectile core without regard to the target is 
to construct an optical potential convolution of a nucleon-ion potential with the target 
density. In the calculations of this work we have determined strong absorption radii by 
this method, using the optical potential of Sec. III. We also take a diffuseness parameter 
a = 0.6/m. 

B. Breakup Amplitude 

In this section we give a more general formulation of the breakup problem which is less 
dependent on the reaction model assumptions of the text. We start with Eq. (1), again 
assuming a straight-line trajectory. Using a Galilean transformation for the initial wave 
function as described in ref.|p, Eq. (1) becomes 

= ^ / dhndz<j)f{r)Vnt{r)e"'f^'4>^{hn - be, h), (19) 

where = ki = —{si — + \mv'^) / {Tiv) and kf^ = k2 = —{si — Sn — \mv'^) / (hv) are the 
2;-components of the neutron momentum in the initial and final state respectively. 
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Introduce now a T-matrix in a mixed representation as 

T(k„, b„, k2) = ^j dz<Pf{h, z)Vnt{h, zy^^\ (20) 



where k„ is the neutron final (measured) momentum. 

If the eikonal approximation for the final wave function is used 



(b,^) =e-^'i^-''"e-^''^^e^^+-''"'^"'^'^'^''', (21) 



the dz integral in Eq.(^) can be calculated by parts when the initial binding energy is 
small and the final neutron energy is close to the neutron incident energy and the neutron 
scattering angle is small such that k2 — Qz — 0. Then 



hv 



J rfze^('=-'?^)^K*(bn,;^)e'^/-«=^"'^""^^'^"*(^'^'^')'^^' = 1 - e-'^(^"), (22) 



where 

X 



(b„) = ^ r e'^''-'^^'^''Vnt{x, y, z')dz' ^ ^V^tih^, 0), (23) 

nv J-oo nv 



where Vnt is a complex potential whose real and imaginary strengths are negative. We 
finally obtain the breakup amplitude in the eikonal form as 

Af^ = J^ rfb„e-*'i^-^"(l-e-*^(''"))0,(b„-be,A;,). (24) 

Similarly a more general T-matrix in a mixed representation can be obtained starting 
from the partial wave form 

T(k„,k^,/C2) = ^^S^„(2/„ + l){Si^ - l)PiS^n ■ ka), (25) 
on the energy shell such that = k^, and defining 

T(k„,b„,A;2) = y^e-^''-^"T(k„,k^,A;2) 

This is unfortunately off shell but we will do the following high energy approximations 
similar to what we have done above in the eikonal case. In the large angular momentum, 

small scattering angle limit, define b' = ^ , we can substitute the partial wave sum 

with an integral and also use Pi„{kn ■ ka) ^ ^o(|k„ - k„|6') = {27i)-^ f dcjye''^'^^-^^^^''^ 
where q_L is the transverse component of k„ . Then 
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-E,„(2/„ + l)(5,„-l) 



1 4:71 



ihv 2m 2ikn 

^ / '^b'„e-''^-^"(^.„ - 1) / dh^ J ^e->'-(^"-^")0,(b„ - be, K) 

= I db„e-^^-^-(l - 5^j0,(b„ - bd, h). (27) 

This formula combined the improved treatment of the interaction by the TC theory 
with the better treatment of the neutron wave function by the eikonal theory. 
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Table I. Energy dependent optical model parameters. aij=0.387fm, r/=1.368fm, a7=0.3fm 
at all energies. 





Vr 


fR 


Ws 


Wy 


(MeV) 


(MeV) 


(fm) 


(MeV) 


(MeV) 


20-40 


38.5 - 0.1455/ 1.447- 


0.005(£/ - 20) 


1.666 + 0.3655/ 


0.3755/ - 7.5 


40-120 






16.226 - 0.1(5/ -40) 


7.5 - 0.02(5/ -40) 


120-180 






8.226 - 0.07(5/ - 120) 


5.9 



Table II. Initial state parameters in Be. Binding energies in MeV, asymptotic normal- 
ization constants Cj in fm~2. 



J, TT 


/ 


j 






1/2+ 





0.5 


3.32 


2.45 


1/2- 


1 


0.5 


3.52 


1.29 



Table III. Breakup cross sections in mb from the initial states 1/2+ and 1/2^ in Be 
on ^Be, at several incident energies. Incident energies inA MeV, Rs in fm. Experimental 
data and shell model spectroscopic factors from ||17|. 



J, TT 




Rs 


O'str 


O'str 




(^diff 


0"-n 


O-n 




SsM 


Stc 


Seik 








TC 


eik 


TC 


eik 


TC 


eik 










1/2+ 


20 


6.5 


35 


31 


45 


46 


80 


77 












40 


6.1 


41 


51.6 


31 


35 


72 


86.6 












60 


5.8 


43.5 


54.6 


24 


31 


67.5 


85.6 












78 


5.6 


43 


54.8 


19 


24.7 


62 


79.5 


32(5) 


0.69 


0.65 


0.51 




100 


5.4 


42 


53 


13 


17.6 


55 


70.6 










1/2- 


20 


6.5 


19 


16.5 


21 


22.7 


40 


39.2 












40 


6.1 


23.4 


29.6 


16.3 


18 


39.7 


47.6 












60 


5.8 


27.2 


32.7 


13 


15.9 


40.2 


48.6 












78 


5.6 


28 


34 


11 


13 


39 


47 


18(3) 


0.58 


0.54 


0.45 




100 


5.4 


28.4 


33.7 


8.4 


9.7 


36.8 


43.4 
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Figure Captions. 

Fig.l Ratio of phase space integrals witli and witliout momentum cutoff, for a (i-wave 
neutron wave function. Tlie effect of the cutoff is to include less than 90%, between 90% 
and 95% and more than 95% of the initial momentum distribution as marked on the 
figure. 

Fig. 2 Top left: Experimental cross sections for neutron scattering on ^Be together with 
the optical model cross sections using the potential of Table I. The full curve is the total 
cross section, the dotdashed curve is the elastic cross section while the dashed curve 
is the reaction cross section. Bottom left: same quantities calculated in the eikonal 
approximation. Top and bottom right: the same quantities calculated with the Varner 
potential. 

Fig. 3 The behavior of the term |1 — S'p and of the transmission coefficient 1 — jS'p 
calculated by the optical model (solid and dashed line respectively) and by the eikonal 
approximation (dotted and dotdashed line). Results are shown as a function of the angular 
momentum /„. In the case of the eikonal calculations the semiclassical relation Z„ + 1/2 = 
kbn was used to make the connection with the neutron impact parameter. Results are 
given at incident energies ranging from 20 to 100 MeV. 

Fig. 4 Total breakup probabilities from the l/2~ state in ^"^Be, as a function of the ion-ion 
impact parameter for diffraction obtained in the TC calculation by the solid line. The 
dashed line is for the stripping. Energies are as in Fig. 3. Dotted and dotdashed lines 
are diffraction and stripping probabilities from the eikonal calculation. 

Fig. 5 The integrand function of the diffraction (a), and stripping (b) term of Eq.(|Il|) 
after integration, full curve, obtained from the realistic bound state wave function and 
the corresponding terms, diamonds, in the sum over partial waves of Eqs.(^,^ in the 
case of the incident energy of 78A MeV. Crosses are the results of a calculation of Eq. 
(|TT]) in which the eikonal phase shifts have been substituted by the optical model phase 
shifts according to Eq. (pT]) . All calculations done at fixed impact parameter be = 5.6 fm 
between the projectile and target. 
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